INVERSE PROBLEMS WITH PARTIAL DATA FOR A 
DIRAC SYSTEM: A CARLEMAN ESTIMATE APPROACH 



MIKKO SALO AND LEO TZOU 

Abstract. We prove that the material parameters in a Dirac system 
with magnetic and electric potentials are uniquely determined by mea- 
surements made on a possibly small subset of the boundary. The proof is 
based on a combination of Carleman estimates for first and second order 
systems, and involves a reduction of the boundary measurements to the 
second order case. For this reduction a certain amount of decoupling 
is required. To effectively make use of the decoupling, the Carleman 
estimates are established for coefficients which may become singular in 
the asymptotic limit. 



1. Introduction 

This article is concerned with the inverse problem of determining un- 
known coefficients in a Dirac system from measurements made on part of 
the boundary. A standard problem of this type is the inverse conductivity 
problem of Calderon [4j, where the purpose is to determine the electrical 
conductivity of a body by making voltage to current measurements on the 
boundary. In mathematical terms, if 7 is a smooth positive function in the 
closure of a bounded domain C R", the boundary measurements are given 
by the Cauchy data set 

C-y = {{u\an,jd^u\dn) ; V • {jVu) = in 17, u£ H\n)}. 

Here uIqq and jd^ulQQ are the voltage and current, respectively, on d^l, 
corresponding to a potential u satisfying the conductivity equation in 17 
{diyU denotes the normal derivative). The inverse problem is to determine 
the conductivity 7 from the knowledge of the Cauchy data set C^. 

The inverse conductivity problem has been well studied, and major results 
include [l], [2^, [27] which prove that C-y determines 7 in various settings. 
Less is known about the partial data problem, where one is given two sets 
ri,r2 CI d^l and the boundary measurements are encoded by the set 

C^''^^ = {(ulr^jdMr^) ; ^ ■ il^u) = in n, u£H\n)}. 

There are two main approaches for proving that 7 is determined by Cy^'^^. 
The first approach, introduced in [3j and [15], uses Carleman estimates with 
boundary terms to control solutions on parts of the boundary. The result in 
|15j is valid in dimensions n > 3 and for small sets r2 (the shape depending 
on the geometry of 50), but assumes that Fi has to be relatively large. 
The second approach [10] is based on reflection arguments and is valid when 
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n > 3 and Ti = T2 and Fi may be a small set, but it is limited to the case 
where dfi \ Fi is part of a hyperplane or a sphere. If n = 2, a result similar 
to [TO] but without the last restriction was recently proved in [12j . 

We are interested in inverse problems with partial data for elliptic linear 
systems. In the case of full data (that is, Fi = F2 = dil), there is an 
extensive literature including uniqueness results for the Maxwell equations 
[M], the Dirac system [H], [2S], and the elasticity system [H]. 
However, it seems that partial data results for systems are more difficult 
to establish. The reflection approach is in principle more straightforward 
to extend to systems, and the recent work [5j gives a partial data result 
analogous to [TO] for the Maxwell equations. As for the Carleman estimate 
approach, there is a fundamental problem since Carleman estimates for first 
order systems, such as the ones in L26j, seem to have boundary terms which 
are not useful in partial data results. 

In this paper, we prove a partial data result analogous to [151 for a Dirac 
system. To our knowledge this is the first such partial data result for a sys- 
tem. The proof is based on Carleman estimates, and it involves a reduction 
to boundary measurements for a second order equation. The corresponding 
boundary term is handled by a Carleman estimate for second order systems, 
designed to take into account the amount of decoupling present in the orig- 
inal equation. In the set where one cannot decouple, we need to use the 
first order structure as well. The Carleman estimates need to be valid for 
coefficients which may blow up in the asymptotic limit, in order to obtain 
sufficiently strong estimates for solutions on the boundary. 

Let us now state the precise problem. We consider the free Dirac operator 
in R^, arising in quantum mechanics and given by the 4x4 matrix 



where D = —iSJ and a = (o"i, cr2) C3) is a vector of Pauli matrices with 



Let C be a bounded simply connected domain with C°° boundary, let 
A € C°°{Q.\'R?) be a vector field (magnetic potential), and let q± be two 
functions in C°°(r2; R) (electric potentials). We will study a boundary value 
problem for the Dirac operator 






Cv = P{D) + V, 



(1.2) 



where the potential V has the form 




(1.3) 




where u± E L'^{Q)'^. By [2ll Section 4], the 



boundary value problem 



Cyu = 
u+ = f 



in fi, 
on do. 
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is well posed if is in the resolvent set of i2y, and then there is a unique 
solution u G H^{n)^ for any / G The boundary measurements 

are given by the Dirichlet-to-Dirichlet map 

Ay : H^/\dn) ^ H^/\dn), / ^ u^ion- 

It is known that the map Ay is preserved under a gauge transformation 
where A is replaced by ^ + Vp where pIqq = 0. Such a transformation does 
not change the magnetic field V x j4, and the inverse problem is to recover 
the quantities V x A and q± from the boundary measurements. 

We are interested in the inverse problem with partial data, where the 
boundary information is the map Ay restricted to a subset T C dQ. More 
generally, we can consider boundary measurements given by the restricted 
Cauchy data set 

Cy = {(n+|an, ii_|r) ; u G H'^{0,)^ is a solution of £yn = in 0,}. 

If is in the resolvent set of Cy, then = {(/, Ay/|r) ; / G H^/^{dnf}. 
Again, the set Cy is preserved when A is replaced by A+Vp where p\qq = 0, 
so the inverse problem is to determine V x ^ and q± from Cy. 

We will prove the following partial data result. Let ch(r2) be the convex 
hull of il, and if xq G define the front face of by 

F{xo) ={xedn; {x- xo) ■ iy{x) < 0}. 

If r C dQ, we write = \ T for the complement in dQ. 

Theorem 1.1. Let Q, C be a bounded simply connected domain with 
connected C°° boundary, let ^1,^2 G C°°{0,;Ji?), and let gi,±,g'2,± G 
C°^(0;R). Let T be any neighborhood of F{xo) in dO,, where xq ^ ch(r2), 
and assume the boundary conditions 

Ai = A2 on d^l, 
Qi,± = Q2,± and d^qi^i = duq2,± on dO., 
qi- 7^ on T^. 

If Cy^ = Cy^, then V x = V x and qi^± = q2,± in 0, 

In the full data case (when T = dQ), the inverse boundary problem for the 
Dirac system and the related fixed frequency inverse scattering problem have 
been considered in [9], [H], [18], [21], [26], [28]. In particular, Theorem ll.il 
for full data was proved in [2T] for smooth coefficients and in [26] for Lipschitz 
continuous coefficients. For T = d^l the boundary conditions (ll.4p - (ll.6p are 
not required, but for partial data results based on Carleman estimates as in 
[3], [15] such conditions are usually needed at least on the inaccessible part 
r*^. By suitable boundary determination results and gauge transformations 
as in |26j . we expect that it would be enough to assume ()1.4p only for the 
tangential components of Ai and A2 on and (jl.Sp only on V^. 

The most interesting condition is ()1.6p , which allows to decouple the Dirac 
system at least on some neighborhood of the inaccessible part T^. This 
decoupling is required for the reduction from boundary measurements for 
Dirac to boundary measurements for a second order system, and also in 



(1.4) 
(1.5) 
(1.6) 
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patching the Carleman estimates for first and second order systems together 
to obtain decay for solutions on part of the boundary. 

Let us outhne the structure of the proof. In Section [2l it is shown that the 
assumption Cy_^ = Cy^ along with (jl.4p ~ (jl.6p implies the integral identity 

/ C/2 (Fi - V2)Ui dx = - [ —U;+duU+ dS (1.7) 

where Ui and U2 are any 4x4 matrix solutions of Cvj Uj = in $7, and further 
U = U1—U2 where U2 is a solution of CV2U2 = in Q with U2^+\dn = f^i,+ |ar2- 
The normal derivative d^Uj^ corresponds to boundary measurements for a 
second order equation. 

The matrices Ui and U2 will be complex geometrical optics solutions to the 
Dirac equation, depending on a small parameter h and having logarithmic 
Carleman weights as phase functions. Such solutions were constructed for 
the Schrodinger equation in [15] and for the Dirac equation in [26]. The 
construction relevant to this paper is presented in Section [3l 

The recovery of coefficients is given in Section [H and proceeds by inserting 
the complex geometrical optics solutions Ui and U2 into (II. 7|) and by letting 
/i ^ 0. With suitable choices, on the left hand side one obtains (nonlinear) 
two-plane transforms of the parameters involved, and microlocal analytic 
methods allow to determine the coefficients. The argument is an analog for 
the Dirac operator of results in [7], and also involves ideas from [21], |26j . 

The remaining issue, and also the main contribution of this paper, is 
the analysis in terms of decay in h of different parts of the boundary term 
in ()1.7p . This is done in Section [5l By a Carleman estimate, we may 
estimate d^U-^. by a second order operator applied to C/+. We will apply an 
/i-dependent decomposition of $7 into a set where q2- is not too small (so 
one can decouple) and where §2,- is small, and the second order operator 
will be chosen accordingly. The coefficients of this operator will typically 
blow up when h becomes very small. 

The second order Carleman estimate is given for a phase function which 
is convexified by a parameter e as in [7] and pTS] , but there is the new feature 
that e needs to depend on /i in a precise manner related to the decomposition 
of f2 to obtain sufficiently strong control of constants in the estimate. In 
the set where q2- is small, we also use a Carleman estimate for the Dirac 
operator to obtain the final bounds. 

More precisely, the proof of Theorem ll.il proceeds in several steps. Noting 
that ()1.7p is an identity for 4x4 matrices, the proof begins by looking at 
the upper right 2x2 blocks in ()1.7|) and by showing that V x yli = V x ^2- 
After a gauge transformation one may assume that Ai = A2, and then from 
the upper left and right 2x2 blocks of (II. 7p one obtains that = q2,-, 
and also = ^2,+ at all points where qi- is nonzero. The coefficients q^ 
would be recovered from the lower right 2x2 block of the integral identity, 
but the estimates for this block in the boundary term seem to be difficult. 
However, at this point one has enough information on the coefficients to go 
back to the Dirac equation and use unique continuation, so that the partial 
data problem can be reduced to the full data problem. Then the result of 
|21j shows that qi^^ = 52,+ everywhere, which ends the proof. 
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Finally, we remark that there is a large literature on Carleman estimates 
and unique continuation, also involving logarithmic weights. We refer to [2], 
[13j . |19j for such results for Dirac operators. Inverse problems for Dirac 
operators in time domain are discussed in |17j . 

2. Integral identity 

The following integral identity will be used to determine the coefficients. 
We write {u\v) = f^v*udx, = {u\u), and (n|ti)s = J^v*udS where u 
and V are vectors or matrices in VL, and S is a subset of dVt. 

Lemma 2.1. Assuming the conditions in Theorem I l.H one has the identity 

((Fl - V2)Ui\u2) = -{-^dyUj^\u2,+)T'^ 
91 - 

for any solutions Uj G ($7)'^ of {P{D) + Vj)uj = in Q, where u = ui — U2 
is a function in {H'^ f] Hq{^}))'^ x H^{i})'^ satisfying duU+\r = 0, and U2 S 
H^{Q)'^ is a solution of {P{D) + V2)u2 = in Q with tt2,+ |an = ""i.+ lan and 
"2 -If = ■"1 -If- 

Proof. Note that the existence of U2 with the stated properties is ensured 
by the condition Cy^ = Cy^ . We first show that 

{{Vi - V2)ui\u2) = i{{a ■ _ - U2-)\u2,+)r'= (2.1) 

Since {P{D)wi\w2) = iwi\P{D)w2) + i(P(i^)u'i 1^2)90 and = V2, we 
have 

((^^1 - V2)UI\U2) = -{P{D)UI\U2) + {m\P{D)u2) 

= i{P{iy)ui\u2)an 

= i{P{v){ui - U2)\u2)dii + i{P{i^)u2\u2)dn- 
Then (j2.ip follows since (ui — U2)+\dn = 0, {ui — n2)-|F = 0, and 

l{P{l^)u2\u2)dn = {U2\P{D)U2) - {P{D)U2\U2) 
= (^2^21^2) - ('U2|V2n2) 

= 0. 

Now u E H^{n)^ with -Au = P{D)[-ViUi + ^2^2) e L'^i^f, and since 
G Hq{Q)'^ we obtain G H^(i^)^ by elliptic regularity. It remains to 
show that 

iqi-{a ■ v)u- = —d^u^ on dQ. (2-2) 
Since ui and U2 are solutions, we have 

a ■ {D + Ai)ui^+ + gi,_ui,_ = 0, 
a ■ {D + A2)U2,+ + q2,-U2- = 0. 

This shows that 

qi-Ui- - q2-U2- = -a ■ Du+ - {a ■ Ai)ui^+ + {a ■ ^2)7x2,+ in fl. 
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Restricting to dO, and using the boundary conditions on the coefficients, and 
writing Du+ = —i{dyUj^)v + (Z)ii_|_)tan on the boundary, where ^tan is the 
tangential component of a vector field A, we obtain 

qi^^U- = —a ■ Du+ = i{a ■ v)duUj^ — a ■ (L'u+)tan on di^. 

Since n+ = on we have (-DM+)tan = on dfl, and (j2.2|) follows upon 
multiplying the last identity hy i{a • i^) . □ 



3. Construction of solutions 

The recovery of coefficients will proceed by inserting complex geometrical 
optics solutions ui and U2 into the identity in Lemma l2.ll These solutions 
depend on a small parameter h > 0, and have the form u = e~P^^m where 
p is a complex phase function and m has an explicit form when /i ^ 0. 

For second order elliptic equations, complex geometrical optics solutions 
go back to [3], [27] in the case where p is a linear function, and they have been 
used extensively in inverse problems for different equations (see the surveys 
[29] , [30] ) . A more general construction was presented in [15] , allowing phase 
functions p = (p + iip where (/? is a so called limiting Carleman weight and 
tp solves a related eikonal equation. See [6] for a characterization of the 
limiting weights. In [I5j . the logarithmic weights (p{x) = log|x — xo\ were 
used to obtain results in the inverse conductivity problem with partial data. 

For the Dirac system considered in this article, a construction of complex 
geometrical optics solutions was given in This construction, special- 

ized to logarithmic Carleman weights, will be reviewed here. Let A, q± be 
coefficients in C°°(0). Instead of 4- vector solutions we will use 4x4 matrix 
solutions U (so that every column of ?7 is a solution) to CyU = in $7, 
having the form 

U = e-P^^{Co + hCi + h^R). (3.1) 

Here /i is a small parameter, p = ip+itp is a complex phase function satisfying 
the eikonal equation (Vp)^ = 0, Co and Ci are smooth matrices with explicit 
form, and ii is a correction term. 

We move to the specific choices of p and Cj, following [7] and [26]. Fix 
a point xq £ \ ch(r2), where ch(0) is the convex hull of $7, and let 
<f{x) = log |x — xo|. We choose 

■ip(x) = distc.2 ( ^,u; 

\\x - xo\ 

where w G S*^ is chosen so that is smooth near Q. Then p = Lp + iip satisfies 
(Vp)2 = near H. 

It will be convenient to make a change of coordinates as in [7]. Choose 
coordinates so that xq = 0, u; = ei, and C {xs > 0}. Write x = (xi,re*^) 
where r > and 6 € (0, vr), and introduce the complex variable z = xi + ir. 
Also write Cr = (0, cos 0, sin 0) and C = + ier- In these coordinates one 
has 



log z, Vp = k, Ap 



z{z 
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As in [26', Section 3], the matrices Co and Ci will be chosen to satisfy 
transport equations involving the Cauchy operator ( ■ D. We will also use a 
function (p G C°°(0) solving 

C • {V(t> + ^) = inn. 

A particular solution cj) is obtained by extending A smoothly into R'^ as a 
compactly supported vector field, and by letting (p = {(^■'V)^^{—(- A) where 
the Cauchy transform is defined by 

{C.V)-^f{x) = ^ [ —^f(x-yiReC-y2lmOdyidy2. 
2vr 7^2 yi + iy2 

Below, we will always understand that (p is this solution. The extension of 
A outside O will not play any role in the final results. 

The following proposition gives the existence and required properties for 
complex geometrical optics solutions. We use the notation introduced above, 
and the notation 

We also write A < B to denote that A < CB where C is a constant which 
does not depend on h. 

Proposition 3.1. Let a e C^iQ) satisfy (C • V)a = in 17. Then for /i > 
sufficiently small, there exists a solution to CyU = in 17 of the form ()3.ip 
where p = log z, 



with C ■ (V0 + A) = in and 

Ci = -{P{D + A)- Qi){r-'/\^^a) + -P{Q)C^ 
I z 

with 11(511154/1,00(^2) < 1. Further, we have 

ll-R||//i(f7) ^ 1- 

Proof. To obtain the H^{Q) estimate for R, in fact we need to compute more 
terms in the asymptotic expansion in terms of h and look for a solution of 
the form 

U = e-p/^{Co + hCi + h'^C2 + h^C-i + h^R^). 

With the choices of smooth matrices Cj given below. Proposition 3.1 in |26j 
implies the existence of such a solution with ||Cj||5yi,oo < 1 and ||-R4||L2(f^) + 

Vi?4||^2(Q) < 1 if /i is small enough. We then obtain the required solution 
(l3TT]) upon taking i? = C2 + /1C3 + hR^. 

The conditions for Cj in [26, Proposition 3.1] are 

Co = P(Vp)Co, MaCo = 0, 
Ci = \{P{D + A)- Qi)Co + PiVp)Ci, MaCi = iHa^wCo, 
C2 = \{P{D + A)- Qi)Ci + P{Vp)C2, MaC2 = iHA,wCu 
C3 = Up{D + A)-Qi)C2. 
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Here Ma and Ha,w are the transport and Schrodinger operators 
Ma = {2Vp ■{D + A) + \Ap)h, 
Ha,w = {D + Afh + 



a ■ iy X A) — q+q-h — o" • Dq+ 

—a ■ Dq^ a ■ (y X A) — q+q-h 



Also, Cj are smooth matrices in solving the transport equations. 
Let (j) and a be as stated. Using the special coordinates, we have 



1 



Ma = -[2C-{D + A) + -]L 



Then Cq = r~^/^e*'^a/4 solves MaCq = in fi, and Cq has the desired 
form. Now one can solve the transport equations C2 and C3 by the Cauchy 
transform for instance, and this shows that also Ci is as required. □ 

Remark. It is possible to perform the above construction of solutions with 
p replaced by —p or p, since these functions also solve the eikonal equation. 
The corresponding forms for the solutions are, respectively. 



U 



U 



1 



h 



-p(C)Ci + oih^) 

z 



-p/h 



h 



LZ 



P{C)r~''^e'^a + -(P(D + A) - Qi)(r~'''e'*a 
i 

+ ^P(c)Ci + oih^ 



where C • (V0 + A) = and C • Va = in and ||C'i||vi/i,oo(q) < 1. 



4. Uniqueness proof 

In this section, we give the proof of Theorem 11.11 modulo the estimates 
for boundary terms which are contained in Section O The following simple 
algebraic identities, valid for a,b £ C^, will be used many times in the 
computations below: 

{a ■ a){a ■ b) + {a ■ b){a ■ a) = 2{a ■ b)l2, {a ■ a)^ = (a • a)l2, 
P(a)P(6) + P{b)P{a) = 2(a • b)h, P{af = {a ■ a)h, 
P{a)Q = QiP{a). 

Since C • C = 0, we also have {a ■ Cf = Q and P{(^f = 0. 

The starting point for the recovery of the coefficients is Lemma [2. 11 which 
implies that 

{{Vi - V2)Ui\U2) = -{—d,U+\U2,+)r^ (4.1) 

where Uj are 4x4 matrix solutions of C-VjUj = in fi, [/ = U\ — U2, and 
U2 solves CV2U2 = in 17 with fZ+jg^ = 0, U-\y = 0. 
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We use Proposition 13.11 or more precisely the remark after it, and choose 
solutions Ui and U2 with 



1 



z 



where ( ■ (Vcpj + Aj) = and ( ■ Vaj = in and where ||-Ri||Hi(n) ^ ^ 

and \\R2\\m{n) ^ h. 

The next result, whose proof is given in the next section, takes care of 
part of the boundary term in (j4.ip . 

Lemma 4.1. The upper right 2x2 block of {-^^dyUj^\U2^+)T'^ is o(l) as 

It is now possible to show that the magnetic field is determined by partial 
boundary measurements. 

Lemma 4.2. V x j4i = V x in Jl. 

Proof. Since P{C)QjP{0 = P{C)P{C){Qj)i = 0, the left hand side of gl|), 
with the above choices for Ui and U2, becomes 

/ U^{Vi-V2)Uidx = - / P{C)P{A,-A2)P{C) ^^dx + Oih). 

The identity P(C)-P(A) = -P(A)P(C) + 2(C • A)l4 implies 

/ U^iVi - V2)Ui dx = -2 P(C)(C • (^1 - ^2)) dx + 0{h). 

Jn Jn z r 

Taking the limit as /i — > in the upper right 2x2 block of (j4.ip . gives 

by Lemma |4. II that 

gi(</.i-</.2)(^ • C)(C • (^1 - A2))aia2Z-'^r-^ dx = 0. 

n 

We choose 01(2:,^) = z'^g{z)b{9) and a2{z,9) = 1, where g{z) is a holomor- 
phic and smooth function in the closure of fig = G C ; (3;i,re*^) G il}, 
and b{9) is any smooth function. Note that C = C(^)- Moving to polar 
coordinates in the x' variables and by varying b{9), we obtain that for all 9 

(a • C) [ e'^^'-^^\Ai - A2) ■ (ei + ier)giz) dzAdz = 0. 
Jng 

Since a ■ is not zero for any 9, it follows that 

gi(</.i-^2)^^^ - A2) ■ (ei + ier)g{z) dz A dz = 0. 

ng 

The last expression is related to a (nonlinear) two-plane transform of 
V X (Ai — A2) over a set of two-planes. We may now apply the arguments 
in [TJ Section 5] (see also [Ml Section 7] , where the last identity is the same 
as formula (40)). One first shows by complex analytic methods that the 
identity remains true with e*^''^!"'^^) ^^id g replaced by 1. It follows that 



/ C-{Ai- A2) dzAdz = 
Jng 
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whenever ^ is in the two-plane spanned by ei and e,.. Varying xq and oj in 
the construction of solutions slightly, this implies that 



[ e • (^1 - A2) ds = o 

Jpnn 



for all two-planes P such that the distance between the tangent space T{P) 
and the point (0, ei) is small. Finally, an argument involving the microlocal 
Helgason support theorem and the microlocal Holmgren theorem shows that 
V X Ai = V X A2 in n. □ 

Since is simply connected and V x ^1 = V x yl2, we see that Ai — A2 = 
Vp for some function p G C°°{fl). Also, by the assumption that Ai = A2 
on do,, we see that p is constant on the connected set dO,. Thus, we can 
assume that p\qq = by substracting a constant. Then Cy^ is preserved 
under the gauge transformation A2 A2 + Vp, and consequently we may 
assume that Ai = A2- We shall write A = Ai = A2 and (j) = (pi = 4>2- 

By Proposition 13.11 there exist solutions Ui and U2 to the equations 
CvjUj = in (j = 1, 2), such that 



1 



+ -{P{D + V(t) + A)- Qi,i){r-^'^ai) - -P{OCi + h^Ri 
I z 

'\p{0r-^''a2 

+ -{P{D -V<P-A)+ Q2,/)(r-^/'a2) + -^2*^(0 + h^R2 

I z 



where (V0 + A) =0, C • Va^ =0, ||Cj ||vt/i,oo < 1, and < 1. 

With these choices for \J\ and U2, we have the following result for the 
boundary term in ()4.ip which will be used in recovering the electric poten- 
tials. Again, the proof is deferred to the next section. 

Lemma 4.3. The upper left and right 2x2 blocks of (^^c?i/C/+|f72,+)r= are 
o(/i) as /i ^ 0. 

From the upper left and right 2x2 blocks of (j4.ip . it turns out that one 
can recover g_ everywhere and qj^ at those points where g_ 7^ 0. 

Lemma 4.4. One has q\ - = ^2,- in 0- Also, = 52,+ at each point of 
Q where q\ - is nonzero. 

Proof. We introduce the notations Q = Qi — Q2 and dj = r~^^'^aj to make 
the formular shorter. Now Vi — V2 = Q, so (|4.ip becomes 

/ \-P{C)d2 + -{P{D -VcP-A) + Q2,l)d2 + -C*2P{0 

Q 



+ h^R2 



-P(C)ai + -{P{D + V0 + A) - Qij)di 
z I 



h 



--P{OCi + h^Ri dx = -{—d^U+\U2,+)r^. (4.2 



z 



1 



91 - 
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Since P{C)QP{C) = P{C)P{C)Qi = 0, the term on the left of dO]) which is 
0(1) with respect to h vanishes. Also, for similar reasons, all terms involving 
fPCOCi and |C'|P(C) and Rj behave like Oih^). Thus we obtain 



P{C)Q{iPiD + V(I) + A)- Qi,/)ai}a2 

- {{PiD -V(P - A) + Q2j)a2}QPiC)ai\ dx + 0{h 

, 1 



^ Jq z 



91, 



9,C/+|C/2,+)r- (4.3) 



Also the terms involving Vcj) + A vanish because 

P{OQP{V(t) + A) + P{V(j) + A)QP{C) 

= [P{OP{V(t> + A)+ P{V4> + A)P{C)]Qi 

= 2[C • (V</) + A)]Qj = 0. 
The expression (|4.3|) becomes 



h 



1 



^ Jn z 



P{C)Q{{P{D) - Qi,/)ai}a2 - {{P{D) + Q2,i)d2}QP{C)ai 



dx 



d,U+\U2,+)r^. (4.4) 



Note that —QQij — Q2Q1 = qh where q = 52, +92,- — qi,+qi,-, so (14. 4p can 
be written as 



h 



1 



« Jn z 



P{C)P{Ddi)d2 - P{Dd2)P{C)di Qidx 



+ 



h 



« Jn z 



-P{C)qdid2 dx + 0{h^ 



91, 



-5,[/+|C/2,+)r- (4.5) 



Now, in the second integral on the left of ()4.5p . the upper left 2x2 block is 
zero. Thus, multiplying (j4.5p by and taking the limit as /i — > in the 
upper left 2x2 block, we obtain from Lemma 14.31 that 



1 



{a • C)(a • Ddi)d2 - {a ■ Dd2){a ■ C)ai (ft,- - 92,-) dx = 0. (4.6) 



In z L 

At this point we make the choices 

oi = r~^/^z6i(0), 02 = r~ 
where hi{6) is a smooth function. Since 



1/2 



dhi 



Vai = -^r-3/2^6ie, + r-^/'^QW + r^^/a^ 

2 at/ 



where = (0, — sin 6, cos 0), we have 



(cr • C)(o" • -Dai)a2 - (d • Dd2){(J ■ C)di 

= -^{a . C)(a • ee)r~'^ - i [(a • C)(^ • e.) - {a ■ e.)(a • C)]r-\. 
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Using the identity 

(a • a)(cr • 6) = (a • b)l2 + ia ■ {a x b), a,b e C^, 
we obtain from (14.611 that 



• - • ee)bi ((71,- - (/2,-)r-2 dx = 0. (4.7) 

Using the condition (jl.Sp . we may extend — g2,± by zero outside 
and therefore we can assume that qi^± — 52, ± £ C^CR^). We write (j4.7p as 

Note that ( and eg only depend on 9. Integrating by parts in 6 and using 
that ^ = iee, we obtain 

I\c7 • O&i ( y^ ^)^"' ^^a^i dr) de = 0. 

Varying bi, it follows that 

(^•0/ ^^'^'''~'^''-\ xi,r,9)r-'dzAdz = 0, 

for ah 9. 

Since cr • C is never zero, we finally get 

for all 0. This implies the vanishing of a Radon transform on certain planes. 
Now varying the point xq in the definition of the direction u; G 5^ in the 
definition of of "0, and varying 9, we obtain from the microfocal Helgason 
and Holmgren theorems as in [7] that 

d{qi- - q2~) -1 n ■ n 
d9 ' 

Thus qi- — q2- is independent of 9. Since qi^^ — q2- G Cc(R^), we obtain 
= (72,- in as required. 
Finally, we return to (14. 5p and now consider the upper right 2x2 block. 
In the first integral on the left this block is zero, so multiplying by and 
letting /i — > in the upper right block gives by Lemma 14.31 that 



/ -(cr • C,)qdia2 dx = 0. 



By a similar argument as above, we obtain that g = 0. Since qi- = q2-, 
this implies ^1,4. = 52,+ at each point where is nonzero. □ 

We have proved that Ai = A2 and qi- = q2- in and that = ^2,+ at 
any point where is nonzero. The next logical step would be to consider 
the lower right 2x2 block of (14. ip to show that = 52,+ everywhere in 
ft. However, the estimates for the boundary term in this case appear to be 
quite difficult. We will choose another route and reduce the remaining step 
to the full data problem, by using unique continuation. 
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Lemma 4.5. Assume the conditions of Theorem ll.il and assume in addition 
that there is some neighborhood W of in J7 such that 

Ai = A2 = A in W, 

Qi,± = Q2,± = q± in W. 

Then Cy^ = C^, that is, the boundary measurements with fuh boundary 
data coincide. 

Proof. Without loss of generahty, we assume that F'^ is connected (if not 
then argue on each connected piece) . By shrinking W if necessary, we may 
assume that also W is connected and q- ^ in W . 

Let (/,(?) be an element of Cy^, so that there is a solution ui G H^{Q,)'^ 
of CviUi = in Q such that ^1^+ = / and ui- = g on dQ. Since Cy^ = Cy^, 
there is a solution U2 € H^{^)^ of £1/3^2 = in satisfying 

ui^+ = U2,+ on dQ, ui- = U2- on T. 

Set u = ui — U2- Then clearly = on dVl and n_ = on L, and by ()2.2p 
we also have d^uj^ = on T. Furthermore, since all coefficients are identical 
in W and q_ 7^ in VF, we have that m_|_ satisfies 

r (-AI2 + 2{A ■ D)h -jz{^- Dq-)a ■ D + Q)u+ =0 in VF, 

\ u+ = duU+ =0 on VF n F, 

where Q is some smooth 2x2 matrix. 

The last system has scalar principal part, and the unique continuation 
principle holds (this can be seen by applying a scalar Carleman estimate 
to both components of «+, for details see [H]). Since W is connected we 
conclude that u+ = in VF, and consequently djjUj^ vanishes on F'^. Since 
g_ 7^ in VF, the relation (12. 2p again implies that ti_ = on all of dVt. 
We have proved that (/, 5) G , showing that Cy'- C C^. The inclusion 
^ao C C^^ is analogous. ' 12 ^ 

We have proved that all the conditions in the preceding lemma hold, so 
we obtain that C^f = C^f . The uniqueness result in [2T] (or [26]) for the 
full data case then implies that qi^+ = ^2,+ in ^- This ends the proof of 
Theorem 11.11 

5. Carleman estimates 

In this section we prove Carleman estimates and establish Lemmas 14.11 
and 14.31 which allow to take care of the boundary term in the identity (14. 1|) . 
This involves an estimate for dyJJj^ on part of the boundary. To explain the 
strategy, we note that any solution u of Cyu = in satisfies 

a -{D + A)u- + g+u+ = 0, 

a -{D + A)u+ + q-u- = 0. 

Then in the set where g_ 7^ 0, the equations decouple and we see that 
satisfies the second order equation 

a-{D + A) [—a ■ {D + A)u+] - q+u+ = 0. (5.1) 
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We will estimate di,u+ on part of the boundary by using a Carleman estimate 
for a second order equation, as in pjj. However, to account for the set 
where q- is small we need to do the analysis very carefully, cutting off the 
coefficients in a suitable /i-dependent way and also letting the convexification 
depend on h. The details are given in the following result. 

We will use below the notation given in the beginning of Section [21 and 
also the sets dQ,± = {x £ dO, ; ±Vip{x) ■ > 0}. Further, we consider the 
semiclassical Sobolev spaces with norm defined by 

\Mh^^, = 11(1 + {hDfr/^u\\L2^^„), u G Cr(R"), s G R. 
In particular we will consider the case s = 1 with the equivalent norm 

scl 

Lemma 5.1. Let A E C°°(H; R^) and q-,q E C~(H), and let be a 
limiting Carleman weight near Vl. Let < a < 1, and let 



A{x,D) 



q{x) 



2A-D-^ia-Dq^)a-D, ^ < ^^b^. 



2A ■ D, otherwise. 



M'^-Dq^)a-A + q, ^ < V|log/i"|, 



9- 

otherwise. 



There exist constants ho, Co, C, where Cq and C are independent of a, such 
that whenever < h < Hq and when 

^ = ^ + ^Y, e{h) = {Co\logh»\)-\ 
one has the estimate 

_^(||e^/S||2 + We^/f^hDyf) - [ du^\e^'^duv\'^ dS 
^[h] Jan. 

< C\\e'^/''h^{-A + A{x, D) + q)vf + Ch^ [ d^e^^'^d^vl'^ dS, 

Jan+ 

for any v E H'^{Q) with v\qq = 0. 

Proof. The proof follows ideas in [15] (see also [26] )■ First let Cq be a 
fixed number, and initially choose ho = 1. Below we will replace /iq by 
smaller constants when needed, and M > 1 will denote a changing constant 
depending only on ip and the coefficients A, g_, and q. 

Write (f = f{ip) where /(A) = A + ^^^i and introduce the conjugated 

operator Po,<^ = ef^^{—l?lS)e~'^l^ = A+iB where A and B are the formally 
self-adjoint operators 

A = {hDf - {Vtff, B = V'fohD + hDoV<f. 

Then, if v is as above, integration by parts gives that 

\\Po,^v\\^ = \\Avf + \\Bvf + ii[A, B]v\v) - 2h\id^^)d^v\d^v)Qn. 
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In terms of symbols one has i[A, B] = /iOp^({a, b}). The hmiting Carleman 
condition imphes, as in yLSj Section 3] and \26\ Lemma 2.1], that 

{a, b}{x, = ^f'i^f\^^\^ + mix)a{x, + Kx, O^x, 

where 

Here we have chosen ho so that max(sup \ 1) < 1/2 for h < ho, which 
ensures that /'((/?) > 1/2. Quantization gives 

4/,2 h 

i[A, B] = ^TTf'{^?Nip\^ + -[moA + Aom + LoB + BoL]+ h^qix) 
e[h) 2 

where g is a smooth function whose norms are uniformly bounded in h. 
Since |V(/7| is positive near il., we have 

- Mh\\v\\ \\Av\\ - MhWvWjji \\Bv\\ - Mh^Wvf - 2h^{{d^ip)d^v\d^v)dn. 

scl 

This used integration by parts and the fact that v\q^ = 0. We obtain 

l|Po,.^'f > + + ^IMP - M'>^ll»ll?,i, 

- 2h^{{d„'f )d^v\d,yv)gQ. 

For the term involving ||f || we note that 

scl 

\\hDvf = {{hD)'^v\v) = {Av\v) + {\V(p\'^v\v) < \\Avf + M\\vf, 
and by the assumptions on /iq we have that 

Therefore, the Carleman estimate becomes 

Next, consider the operator 

= h^e^/''{-A + A{x, D) + q)e-^'^ 
= Po,(p + hA{x, hD + iV(^) + h^q. 

We have 
Me{h) 

4|| ~„.||2 



\v\\^i^ - 2h-\{du(f)duv\duv)dn 



scl 



< 4\\P^v\\^ + Ah'\\A{x, hD + iV(p)v\\' + 4h^\\qv\\\ 
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If ho is chosen so that |log(/i")| > 1 for h < ho, then by the definition of 
A{x, D) and q{x) it holds that 



\A{x, hD + iV<:p)v\\ < Mv^|log(/i")|||u||^i 



scl 



\\qv\\ < My^|log(/i")|||?;||. 

Thus 

-Jf—\\v\\l, -2h\{d,^)d,v\d,v)9n<M\\P^vf + Mh^\log{h'')\\\v\\l, . 

IVlEytl) scl scl 

At this point we choose Co so that Cq > 2M^, which implies 

M/i2|iog(/i")| < 1— . 

' ^' - 2Me{h) 

With this choice, we arrive at the Carleman estimate 

^^—\\v\r 

scl 



v\\l, - 2h^{{d^^)d,.v\duv)dn < M\\P^vf. 



Me{h) 

Finally, replacing v by e'^^^v gives that 
Me{h) ' 



\e^l\f + \\hD{e^/^v)f) - 2h\{d,^)e^'''dMe^'''d,v)a^ 

< M\\e^'^h^{-/\ + D) + q) 



|2 



Since dyip = f'{ip)dyLp, the result follows. □ 

Remark. The motivation for the choice of e{K) comes from the fact that 
^v/h <; j^-Ca^ip/h -^^j^gj-g (J jg independent of a. Here the factor /i"*^" can 
be controlled by positive powers of /i if a is chosen small enough. 

It will be essential to use the convexified weight (p instead of since we 
will need the stronger constants obtained from convexification to carry out 
the estimates for boundary terms. Next we give a Carleman estimate for 
the Dirac operator, which will also be required for controlling the boundary 
terms. 

Lemma 5.2. Let A £ C°°(Tl; R^) and q £ C°°{Tl), and let if he a limiting 
Carleman weight near 0,. Suppose a, Cq, e{h), and (p are as in Lemma l5. II 
There exist C, Hq > 0, with C independent of a, such that for < /i < /io 
one has 

\\e^/''uf < Ce{h)\\e^/^{a ■ {D + A) + ql2)uf, u E Hl{nf . 

Proof. We follow the argument in [265 Lemma 2.2], where more details are 
given. The Carleman estimate in Lemma 1 5 . II implies that 

/l2 



^1 < C\\e^/^{-h^M2)e^'^u\\l2, 



E{hy 

for all u G C^(i7)^. It is possible to shift the estimate to a lower Sobolev 
index and prove that for h small one has 

/l2 



£{fl) scl 
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Now write e^/^{-h'^ /^l2)e^/^ = P| where P^{hD) = e^/^(a ■ hD)^^'^. 
Since {hD)~^ P^{hD) is an operator of order 0, we obtain 



e{h) 



u\\l2 < C\\e^'^{a-hD)e-^'^u\\l2. 



If h is small (so l/£{h) is large), we may replace a-hD by a ■ {hD + hA) + hql2 
in the last inequality. This shows the desired estimate for u G C^(r2)^, and 
the result is then valid for u G Hq{Q)'^ by approximation. □ 

We may now give the proof of Lemma l4.ll which provides an estimate for 
the part of the boundary term required for determining the magnetic field. 

Proof of Lemma Recall that Uj are solutions of Cy.Uj = in fi, where 



p/h 



Ui,+ = e 



TJ* — p-p/h 



1 



-r 



-1/2 i4, 



z 



2x4 



'^■^/4x2 V^2/4x2 



with ||-Ri||/fi(n) ^ h, \\R2\\m{n) ^ h. Also, U = Ui — U2 where U2 solves 
CV2U2 = in with U2,+ \dn = ^i,+ |sn- Thus ?7+|ar2 = 0, and we have 
U-\t = duU+\T = and also U+ G {H"^ n HI{Vl)Y'^^ by Lemma EH 

Denote by J the upper rig ht 2 X 2 block of {-^duU+\U2,+)r-- Writing 

W = Wi — W2 where Wi and W2 are the right 4x2 blocks of Ui and f72, 
respectively, we have 

j= [ e-p/'^Ro—d,W,dS. 



Qi- 

Since qi- 7^ on F'^, we get (using the Frobenius norm on matrices) that 
Jf<c( I \\R2fdj!\( I \\e-^'^dyW+fdS 



Write —ip = — (/? + for the convexified weight corresponding to —(p, 

as in Lemma EH We note that e~^f^^ = me~'^l^ where < m < 1, and also 



the estimate ||i?2||L2(r<=) ^ C'||i?2||//i(n) < h which follows from the trace 
theorem. These facts yield 

ll^f </i'l|e-^/"9,W^+lli2(rc). (5.2) 

At this point we wish to use the Carleman estimate of Lemma l5.ll This 
allows to estimate d^W^ by a second order operator applied to W+. Since 
= — W2,+ where Wi is a solution with explicit form which can 

be estimated, we want to choose the operator so that it will make terms 
involving W2,+ vanish. Note that when ^2,- 7^ 0, W2,+ solves (columnwise) 
the equation 

a-{D + A2) (^'^ • (D + ^2)^-2,+^ - q2,+W2,+ = 0. 
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In the set where q2- ^ this may be rewritten as 

[-A/2 + 2(^2 • D)l2 - —{a ■ Dq2,-)a ■ {D + A2) + q2h]W2,+ = 0, (5.3) 

where q2 = A2 • A2 + D ■ A2 — g2,+92,-- Since VF+I^q = 0, the estimate in 
LemmaEHl apphed to —if and A2, q2,-, and q2, shows that 

-^{\\e-'^/''W+f + \\e-'^/''hDW+f) + h^ [ OMe'^^^^d^W+f dS 
^(h) Jdn+ 

<\\e-^/''h\-A + A2{x,D) + q2)W+f-h'^ [ OMe'^^^d^W+l^ dS. 

Recah that (9,yTy+|r = 0. Now F is a neighborhood of the front face F{xo), 
but we have F{xq) = d^l- since ip was the logarithmic weight. This shows 
that the last boundary integral vanishes. Since d^f > on F*^, the Carleman 
estimate can be written as 

[ \e-'^/''d^W+\^ dS + ^{\\e-^/''W+f + \\e-^/^hDW+f) 

< \\e-^/hh\-A + A2{x, D) + q2)W+f. 
Going back to (|5.2p . we have arrived at 

\\jf + -^{\\e-~*'^W+f + \\e-^/^hDW+f) 
e{h) 

< /i3||e-'^/^(-A + A2{x,D) + q2)W+\\\ (5.4) 



Let = {x G ; — — < y^|log be a subset of Q where (72,- is 
bounded away from zero, with an /i-dependent bound. The proof will then 
be completed by establishing the following two estimates: 

h^e-^/\-A + A2{x,D)+q2)W+\\l,^g^^=o{l), (5.5) 

h^e-^/h{-A + A2{x, D) + g2)Vf^+|li2(^,5,^) 

< 0(1) + hiWe-^/'^W+f + \\e-^l''hDW+f ). (5.6) 

Proof of (1531). By we have (-A + A2{x,D) + q2)W+ = (-A + 

A2{x,D) + (72)1^1,+ in Sh- To obtain sufficient decay in h, we will need 
to convert — AWi^_|_ into first order derivatives of Wi,± by noting that Wi 
solves 

a-{D + Ai)Wi^+ + = 0. 

Then applying a ■ {D + Ai) implies 

-AWi^+ = -2Ai ■ DWi^+ - qi^-a ■ DWi^- + Qi,+Wi^+ + 

for some smooth matrices Qi,±. Thus, the most significant terms in the ex- 
pression {—A+A2{x^D)+q2)Wi^-\-, regarding growth in /i, are A2{x,D)Wi^+ 
and -2Ai ■ DWi^+ and -qi-a ■ DWi-. Since Wi,+ = eP'^M and Wi- = 
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eP/^R where \\M\\hh^^) < 1, \\R\\m(n) ^ h, we have 

h^e-^/\-A + A2{x, D) + q2)W+\\l,^Sn) 

< h\\e-^/'^eP/'A2ix, hD - iVp)M\\l,^s,) + h\\e-^/^eP/^ M\\l,^,^^ 

+ h\\e-^l^ePl^R\\l,^s,) ^ /i'-^"(|log/i°| + 1) 

since e""^/^ < where C is independent of a. Choosing a > so 

smah that 1 — Ca > 0, this goes to zero as /i — > 0. 

Proof of (15. 6p . In J7 \ the coefficient q2- is close to zero, and we will 
use that W2 solves the equations 

a-{D + A2)W2,+ + q2,-W2,- = 0, 

a-{D + A2)W2,- + q2,+W2,+ = 0, 

which implies that 

(-A + 2A2-D + q2)W2,+ + {a ■ Dq2,-)W2,- = 0. 
By the definition of A2{x, D) and q2{x), we have on 17 \ S/i 

(-A + 2i2(x, D) + q2)W+ = (-A + 2A2-D + 92)^^1,+ 

- {a ■ Dq2,-)W^ + {a ■ Dq2,-)Wi^^. (5.7) 

Note that we have written VF2,- in terms of W- and For the first 

term on the right hand side of (j5.7p . a similar argument as in the proof of 
(j5.5p implies 

h^e-'^/\-A + 2A2-D + 92)^1,+ ||i2(n) < = o(l) 

when a is small enough. Since Wi^^ has explicit form, the third term satisfies 
h^e-^/\a . Dg2,-)^i,-|li2(^) < h'-'^''^ = o(l). 
To prove (15. 6p . it remains to show that 

h^e-^/\a . Dq2,-)W^\\l.^^^S,) 

< 0(1) + h{\\e-'^/^W+f + \\e-'^/''hDW+f). (5.8) 

To this end we will apply the Carleman estimate for a Dirac operator given 
in Lemma 15.21 This will allow to estimate W- by o" • (D + A2)W-, which 
again may be broken into terms involving the explicit solutions Wi^± and 
the term W+ which is admissible. 

However, the Carleman estimate only applies to functions vanishing on 
the boundary. One has W-\r = since Ui- = U2,- on T, but W- could be 
nonzero on T^. Here we are saved by the fact that the estimate is over the 
set Q Sfi which has to be a positive distance away from if h is small, 
by the assumption that ^2,- 7^ on T^. Thus, let V be a neighborhood of 
r*^ in which 52,- / 0, choose xo G C^{V) with xo = 1 near V^, and let 
X = 1 - Xo- Then 

h^e-^/^a . Z)g2,-)TV-Ili.(n.5.) ^ h^e-^/\W^\\l,^^y 
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Applying Lemma [Ol to x^- £ -ffo(f2)^^^ gives 

\\e-^"'xW-\?mn) < <me~^'\'J ■ {D + AMxW-)\\l.^^^ 
< eW||e-^/^(a • (D + A^m^Wl^^^^ + eih)\\e-^^\a ■ Dx)W^\\l,^^y 
We write 

a-{D + A2)W- = a-{D + - a ■ {D + ^2)W^2,- 

+ a ■ {A2 - Ai)Wi^_ 

= -q2,+W+ + {q2,+ - qi,+)Wi^+ + a ■ {A2 - 

Thus 

\\e-'^/^{a-iD + A2))W.\\ < \\e~'^/''W+\\ + \\e-'^/^Wi,±\\ 

< \\e-'^/^W+\\ + h^^'^ 

by the exphcit form of Wi. Finahy, since q2- / on the support of a ■ Dx, 
we have in this set 

W- = Wi- + —a ■ {D + A2)W2 + 

92- 

= - —a ■ (D + A2)W+ + —a ■ {D + ^2)^^!,+ 

92- 92- 

= Wi-- —a ■ (D + A2)W+ + —a ■ (A2 - Ai)Wi + - _. 

92- 92- ' 92- 

Consequently 

Combining these estimates gives 

h'\\e-^/\a-Dq2,^)W^\\l.^^^s,) 

< h^-^"e{h) + /ie(/i)(||e-^/'^VF+f + \\e~'^/^hDW+f). 
This shows (I5.8p if a is chosen small enough. The proof is complete. □ 



Next we will prove Lemma 14. 3|, which is used in recovering the electric 
potentials. The stronger decay of suitable blocks in the boundary integral 
{o{h) instead of o(l) as in Lemma |4.1|) is due to the fact that Ai = A2 = A. 
Otherwise, the proof will be mostly parallel to that of Lemma l4.ll 



Proof of Lemma \4-3\ The solutions Ui and U2 have the form 



Z 



ri -i/c^ / \ , /R 



-r ''^a2{ + , 

^/4x2 V-"2/ 4x2 



-Z 

with \\Rj\\m(n) < h, \\R'j\\m{n) < h. 

We denote the upper left 2 x 2 K 
show that J = o{h). The argument for the upper right block is analogous. 



We denote the upper left 2x2 block of {-^dM+\U^)r- by J, and will 
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If Wi, W2 are the left 4x2 blocks of Ui and U2, respectively, and if Vl^ = 
W1-W2, then 

J = [ e-p/^e-''^R2—d^W+dS. 
Jr'= Qi- 

Repeating the argument in the proof of Lemma l4.H we obtain the estimate 

+ -^{\\e"^^^W+f + \\e-'^/''hDW+f) 
e{h) 

< h^e-'^/\-A + A2{x, D) + q2)W+f. 

Here A2 and q2 are the coefficients in Lemma 15.11 for A2 = A, q2,-, and 
q2 = A- A + D ■ A - q2^+q2-- Let Sh = {x £ il.; \-^\ < ^/]ogh^} as before. 
Then the desired conclusion J = o{h) will be a consequence of the following 
two estimates: 

h^e-*l^{-A + A2{x,D)+q2)W+\\l,^s^^=o{h''), (5.9) 

h^e-^/\-A + i2(x, D) + q2)W+\\l,^^^s,) 
< o{h^) + h{\\e-^/''W+f + \\e-^/''hDW+f). (5.10) 
Proof of (j5.9p . In Sh one has 



{-A+A2{x,D)+q2)W+ = {-A + 2A-D —{a-Dq2 -)a-{D+A)+q2)Wi +. 

Using that Wi is a solution of the Dirac system with Ai = A, we obtain 
after some computations that 

(-A + 2A • D)Wi,+ = Mi,+VFi,+ + Mi,_M/i,_ 

in fi, where Mi^± are smooth matrices in 0, which are uniformly bounded 
with respect to h. Using the Dirac equation again, we have in Sh 

i-A + A2{x, D) + q2)W+ = Mi,+Wi,+ + Mi,-iyi,_, 



where Mi^± are matrices in Sh satisfying UMi^-tH^oo^^^) < ^/jlog By the 
explicit form for Wi,±, one has 

h^\\e-^l\-A + A2{x,D) + q2)W+\\l,^s,)^h^~''''\\ogh^\. 

This proves (j5.9p if a is chosen small enough. 

Proof of ()5.10p . In J7 \ S"/! we obtain the identity (j5.7p where A2 = A: 

(-A + 2i2(x, D) + q2)W+ = {-A + 2A-D + q2)Wi^+ 

- {a ■ Dq2,-)W- + {a ■ Dq2,-)Wi,-. 

As in the proof of ()5.9p . it follows that 

h^e-^/\-A + 2A-D + g2)W^i,+ |li2(f,) < /i^^"^" = o{h') 

for a small. Also, clearly 

h^e-^/\a ■ DQ2,-)W^i,-|li2(f,) < /i^-^" = o{h^). 
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Using the cutoff x the Carleman estimate of Lemma 15.21 in the same 
way as when proving ()5.6p . we have 

We note that 

a-{D + A)W- = -q2,+W+ + (^2,+ - m n, 

and 

W- = Wi-- —a 

92 - 

This imphes that 

< h^-^'^e{h) + he{h){\\e-^/^W+f + \\e-^/^hDW+f). 
Now h^~^"e{h) = o{h'^) for a smah, so we have proved (IS.lOp . □ 
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